Successive magnetic phase transitions and multiferroicity in Spin-1 triangular lattice 

antiferromagnet Ba3NiNb209 
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We report the magnetic and electric properties of Ba3NiNb209, which is a quasi-two-dimensional spin- 
1 triangular lattice antiferromagnet (TLAF) with trigonal structure. At low T and with increasing magnetic 
field, the system evolves from a 120 degree magnetic ordering phase (A phase) to an up-up-down (it ltd) phase 
(B phase) with a change of slope at 1/3 of the saturation magnetization, and then to an "oblique" phase (C 
phase). Accordingly, the ferroelectricity switches on at each phase boundary with appearance of spontaneous 
polarization. Therefore, Ba3NiNb209 is a unique TLAF exhibiting both uud phase and multiferroicity. 

PACS numbers: 75.40.Cx, 75.1().Jm, 77.22.Ej, 75.85.+t 
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The two dimensional (2D) triangular lattice antiferromag- 
net (TLAF) is one of the simplest possible geometrically frus- 
trated spin systems When the value of the spin is small, 
5 = 1/2 or 1, in addition to frustration there are strong quan- 
tum spin fluctuations at low temperatures that lead to exotic 
ground states iQ]. One celebrated example is the unusual spin 
up-up-down (uud) ground state. Theories have predicted that 
quantum fluctuations should stabilize the uud state in spin- 
1/2 TLAFs isl-IHt]. In this collinear state, the low temperature 
magnetization is constant over a finite range of the magnetic 
field and equal to 1/3 of the saturation magnetization M^. Ex- 
perimentally, Cs2CuBr4 is a rare example of spin- 1/2 TLAF, 
in which the uud state occurs Recently, the 1/3 Ms 

magnetization plateau has been observed for Ba3CoSb20g, 
another TLAF with effective spin- 1/2 ([§]. The only reported 
spin-1 TLAF exhibiting 1/3 is Ba3NiSb209 0\. Search- 
ing for new TLAFs with magnetization plateaus is still a very 
active topic in condensed matter physics. 

Another interesting property of TLAFs is multiferroicit y, in 
which the magnetic and electric degrees of freedom ifiol ll ill 
are strongly coupled. The multiferroicity has been observed in 
TLAFs with a 120 degree magnetic ordering phase for classi- 
cal spins, such as ACr02 (Cr^+ S = 3/2, A = Ag and Cu) 
and RbFe(M04)2 (Fe3+ S = 5/2) Q. The study on 
RbFe(M04)2 also shows that the ferroelectricity disappears in 
the uud phase, which is stabilized by the thermal fluctuations. 
The theory for multifeiTocity in TLAFs is still under develop- 
ment 11511 , since the well established theories, such as the spin 



current model or the inverse Dzyaloshinshkii-Moriya model 
do not predict ferroelectricity in TLAFs lIlwlSll . However, 
it has been predicted that a broad range of trigonal materials 
with the 120 degree spin structure could be multiferroic, since 
both ACr02 and RbFe(M04)2 have trigonal structure. 

So far, no other TLAFs with trigonal structure and 120 de- 
gree spin structure were found to be multiferroic to further 
confirm this prediction. Also, no TLAFs with quantum spins 
{S = 1/2 or 1) have been reported to show multiferroicity. 
Moreover, the coexistence of the uud phase and multiferroic- 



ity has not been reported in any TLAF. Here, we report de- 
tailed magnetic and electric property studies on new spin-1 
TLAF with trigonal structure, Ba3NiNb209 (Ni'^+ has S = 1). 
We observed successive magnetic phase transitions at 1/3 and 
VS/ 3 of the saturation magnetization, which are also accom- 
panied by spontaneous reversible polarization and dielectric 
anomalies. Although a distinctive magnetization plateau was 
not observed due to the polycrystallinity and the finite temper- 
ature, our results strongly suggest that the title compound is a 
rare example of TLAF with coexistence of collinear magnetic 
structure and multiferroicity. 

Polycrystalline Ba3NiNb209 samples were prepared by 
solid state reaction. Appropriate mixtures of BaCOs, NiO 
and Nb205 were ground together, pressed into pellets, and 
then calcined in air at 1230 °C for 24 hours. The room tem- 
perature powder x-ray diffraction (XRD) pattern [Fig. 1(a)] of 
the as-prepared sample can be indexed as a hexagonal P-3ml 
structure with a = b = 5.7550(5) A and c = 7.0656(2) A as 
reported previously lfl9ll . which belongs to the trigonal space 
group. The 1:2 ordering of the B site ions (-Ni-Nb-Nb-) in 
this structure has been confirmed by two evidences; (i) the 
appearance of low intensity Bragg peaks in the XRD pattern 
related to the ordering structure 12011 . such as the (112) and 
(103) peaks [inset of Fig. 1(a)] and (ii) the high angle annular 
dark field (HAADF) Z-contrast image from scanning trans- 
mission electron microscopy (STEM) [Figs. l(c and d)]. The 
intensity of the Nb and Ni atom columns in the image clearly 
shows the -Ni-Nb-Nb- ordering viewed down the [010] direc- 
tion due to their atomic number Z difference 12 ll 12211 . No de- 
fects are seen on the Ni or Nb sites. We noticed that with the 
annealing temperature above 1230 °C, the sample will turn to 
be Ba(NiQ 33Nbo. 67)03 with B-site disorder and Pm-3m cu- 
bic structure. With the P-3ml structure, Ba3NiNb209 can be 
represented as a framework consisting of corner-sharing NiOe 
and NbOe octahedra [Fig. 1(e)]. The Ni ions occupy the lb 
Wyckoff sites, and this site forms the triangular lattice in the 
ab plane [Fig. 1(f)]. Therefore, the structure can be regarded 
as a pseudo-two-dimensional triangular magnet, i.e., the Ni 




FIG. 1 : (Color online) (a) The room temperature XRD pattern (Cu 
Kai A = 1.54059 A) for Ba3NiNb20g (plus marks). The solid curve 
is the best fit from the Rietveld refinement using FullProf. Inset: the 
(112) and (103) peaks from the XRD pattern, (b) The neutron diffrac- 
tion pattern measured at 2 K showing several magnetic Bragg reflec- 
tions, (c) Atomic resolution STEM HAADF Z-contrast image look- 
ing down [010]. Ba, Nb, and Ni atoms show brightest, brighter, and 
weak contrast, respectively, (d) Line profile of the atomic columns 
intensity indicated by the box in (c) showing the ordering sequence 
of -Ni-Nb-Nb-.(e) Schematic crystal structures for Ba3NiNb209: the 
green octahedra represent Ni sites and the blue octahedra represent 
Nb sites, (f) Zero field magnetic structure below Tni composed of 
Ni^+ ions and the arrows representing their spins. 



magnetic triangular lattices are magnetically separated by the 
two non-magnetic Nb layers. 

The susceptibility x measured at = 0.1 T [Fig. 2(a)] 
shows a peak around 4.9 K related to the long range mag- 
netic ordering. The linear Curie-Weiss fit of \lx at high tem- 
perature (not shown here) gives a 6'cw = -16.4 K, showing 
the antiferromagnetic nature of the exchange interactions and 
an effective moment /ieff = 3.15 /ifi/Ni, which corresponds 
to a Lande g-factor g = 2.Ti and is typical for Ni^+ ions 



1 2311 . Assuming that each Ni^+ spin only interacts via ex- 
change J with the z = 6 nearest neighbor spins on the trian- 
gular lattice, the Heisenberg Hamiltonian is J'Yli^i 
In mean-field theory, this leads to a Curie-Weiss law with 
6'cw = -zJS(SH-l))/3kB. For BasNiNbsOg, 5" = 1 and J/kg 
= -l/(46'cw) = 4. 1 K. The neutron diffraction pattern measured 
at 2 K [Fig. 1(b)] shows magnetic Bragg peaks at Q = {n\ + 
1/3, 1/3, 7i2 + 1/2) {rii. integer). Accordingly, its refinement 
(not shown here) reveals a magnetic structure with collinear 
AFM spins between the nearest neighbor layers and the 120 
degree AFM ordering within the triangular lattice with an or- 
dered moment of 1.8 /ie for each Ni^+ ion. This refinement 



FIG. 2: (Color online) (a) Temperature dependence of x for 
BasNiNbaOg at different fields. Inset: dx/dT vs. Tfoi H = 7 and 
9 T. (b) Temperature dependence of Cp/T at 0, 3, 5, 7, 8, and 9 T. 
Inset: Expansion of the 9 T data, (c) Temperature dependence of e' 
at 0, 4, 8, 10, 12, 14, 16 and 18 T, and (d) magnetic field dependence 
of e' at 0.5, 0.9, 1.6, 2, 2.6 and 3.1 K. Inset: Expansion around H = 
8 T for 0.5 and 1.6 K (data is offset for clarity). The vertical lines in 
(a) and (b) indicate Tni (higher T) and (lower T). 



also shows a B sites ordered P-3ml lattice structure at 2 K, 
which means there is no obvious structural distortion below 
4.9 K. 

For H > 5 T, the peak of the susceptibility x becomes 
broader [Fig. 2(a)] and the related dx/dT vs. T curves display 
two peaks [inset of Fig. 2(a)]. We define Tni and by the 
locations of the higher and lower temperature peaks, respec- 
tively. Below 9 T, both Tni and Tm, decrease with increasing 
field. The zero field specific heat data [Fig. 2(b)] shows a sin- 
gle peak around Tni . For H > 5T, this peak becomes broader 
and could be regarded as two peaks as indicated by lines in 
Fig. 2(b) for 9 T. Their locations are consistent with those of 
Tni and Tn2 defined from dx/dT. 

Samples with rectangular shape with typical dimensions of 
5.0x4.0x0.15 mm-^ were cut from the pressed pellet to mea- 
sure the dielectric constant and the pyroelectric current using 
the techniques found elsewhere |24tl . The real part of the di- 
electric constant (e') measured at zero field [Fig. 2(c)] shows 
a sharp peak at Tni. With increasing field, this peak shifts to 
lower temperatures and becomes broader For i7 > 8 T, a sec- 
ond peak well separated from the first peak appears at lower 
temperatures. The field dependence of e' measured at 0.5 K 
shows two peaks around 8 T and 15 T, respectively [Fig. 2(d)]. 
With increasing temperature, the peak around 8 T does not 
move significantly while the peak around 15 T moves to lower 
fields. The two peaks tend to merge at about 3 K. 

The pyroelectric current was measured on warming after 
poling the crystal in an electric field while cooling down from 
above Tni with magnetic field (FC) or without magnetic field 
(ZFC). The spontaneous polarization was obtained by inte- 
grating the pyroelectric current with respect to time. At zero 
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FIG. 3: (Color online) Pyroelectric current (Ip) (a) and polarization 
(b) as a function of temperature at different magnetic fields under 
positive electric poling field and FC conditions. The inset of (b) 
shows polarization vs. magnetic field at 2.3 K and 3.8 K under ZFC 
condition. Ip (c) and polarization (d) at H = and 10 T under pos- 
itive (solid symbols) and negative (open symbols) poling fields and 
ZFC and FC conditions. 



magnetic field, the pyroelectric current, and accordingly the 
spontaneous electric polarization, begin to develop below Tni 
[Fig. 3(a)]. With opposite poling electric field, the pyroelec- 
tric current direction can be reversed [Fig. 3(c)]. This indi- 
cates the ferroelectric (FE) nature of the ground state, which 
is clearly coupled with the 120 degree magnetic ordering. The 
obtained polarization is around 1 . 1 jiClw? at T. The polar- 
ization overall decreases and shifts to lower temperatures with 
increasing magnetic field [Fig. 3(b)], which is consistent with 
the behavior observed for the higher temperature peak in e' . 
For the polarization emerging at Tni, there is no difference 
between the ZFC and FC measurements. 

For i7 > 8 T another peak of the pyroelectric current (or 
another step-like increase of the polarization) is found at the 
same temperature where the lower temperature peak of e' is 
observed. Note that this low temperature polarization increase 
is observed only when the sample is cooled under the presence 
of both, magnetic and electric, fields, i.e., after magnetoelec- 
tric (ME) annealing. For example, the lower temperature peak 
for the pyroelectric current disappears for the ZFC measure- 
ment, as seen in Fig. 3(c) for the 10 T data. For some ME sys- 
tems, the reversible pyroelectric peak is also observed under 
ME annealing without invoking a ferroelectric transition as a 
consequence of the combined effects of magnetic domains and 
the ME effect |24-2^. Therefore, the additional increase of 
the polarization and the low temperature e' anomaly are more 
likely to arise from a magnetoelectric mechanism, rather than 
from another FE transition with a different order parameter 

Combining all the values for the transition temperatures and 
critical magnetic fields observed from the various techniques 
described above, the low temperature H-T phase diagram for 
Ba3NiNb209 shown in Fig. 4(c) is obtained. Below 4.9 K and 



FIG. 4: (Color online) (a) DC magnetization {M) and its deriva- 
tive (AMIAH) for T = 1.4 K. (b) AC susceptibility as a function of 
magnetic field H for four temperatures. The asterisks and the tri- 
angles indicate A-C and C-D phase boundary respectively, (c) H-T 
phase diagram for Ba3NiNb209. Different symbols denote phase 
boundaries obtained from different techniques. □: e' ff-sweep, •: 
e' T-sweep, o: polarization, <: Cp, A: AC susceptibility ff-sweep, 
A: AC susceptibility T-sweep, ©: DC susceptibility Tni, DC 
susceptibility Tm, *'■ DC magnetization. The arrows indicate spin 
structures.The dashed (solid) vertical lines in (a) and (b) indicate the 
transition fields to uud (oblique) phases. The dashed lines in (c) are 
guides to the eyes. 



with increasing magnetic field, at least three magnetic phases 
(labeled as A, B, and C in Fig. 4(c)) are consistently observed. 
In the A phase, the spins form the 120 degree coplanar AFM 
ordering at zero field. Under small external magnetic fields, 
the order in the A phase is expected to gradually deviate from 
the 120 deg structure while maintaining a similar degree of 
continuous degeneracy 12711 . At higher fields, as the theories 
predict jsl-Sj, for a AFM triangular lattice, the quantum fluc- 
tuations for small spins can lift the classical degeneracy in the 
A phase to select the B phase. We tentatively assign the uud 
state with a 1/3 M» magnetization plateau stabilized in a finite 
magnetic field range to the B phase [3-5S] . 

In order to examine the possibility of the magnetization 
plateau in the B phase, the DC magnetization and AC sus- 
ceptibility (x') on were measured. As seen in Fig. 4(a), the 
DC magnetization saturates around 23 T with Mg = 2.2 /ie- 
A slope change of M vs. H around 8 T is observed, at which 
the magnetization is 0.77 ji-g w 1/3 Ms- Experimentally, the 
AM I AH curve shows a distinctive feature characterized by a 
broad valley in the magnetization plateau region and a sharp 
peak and a broad shoulder respectively on entering and exit- 
ing the putative plateau region Tsi 191 12811 . In Ba3NiNb209, the 
most distinctive feature (valley) is not apparent in the AM I AH 
curve as a consequence of the lack of the plateau-like mag- 
netization. On the other hand, the x' measurement, being a 
technique to probe AM/AH directly, revealed all the features 
as shown in Fig. 4(b): a dominant peak at 8 T followed by a 
broad valley at 9 T and a broad shoulder like peak at around 
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14 T. However, the x' at the valley is not the smallest as the 
overall x' increases with field. There are two possible reasons 
why the magnetization plateau was not clearly observed: (a) 
the polycrystalline nature of our sample and (b) that even in a 
single crystal the plateau is manifested only in certain direc- 
tions of the field with respect to the crystalline axes. 

The 14 T feature is also seen in AM I AH as a broad 
peak, where the magnetization is approximately \/3/3 Ms 
[Fig. 4(a)]. This suggests the 14 T feature is due to the transi- 
tion from the uud to a coplanar 2: 1 canted "oblique" phase (C 
phase) with one spins in the A phase rotated to be parallel with 
another spin, which gives V3/3 Ms- We emphasize here that 
our spin structure model in the oblique phase is based on the 
magnetization data and more rigorous experiments like neu- 
tron diffraction should be performed to confirm the structure. 
In general, the spin structure in the oblique phase, a generic 
term for a spin structure with two parallel spins tilted with re- 
spect to the third, is unique to each TLAF system depending 
on the magnitude of the spin and hence on the order of the 
fluctuations. It is noteworthy that our suggested 2:1 canted 
spin structure is the same as that of RbFe(M04)2 under 10 T 
in-plane magnetic field, which was confirmed by neutron ex- 
periments (see Fig. l(d)ofRef. lll4ll ). At higher temperatures, 
the x' data also show features (indicated as asterisks and tri- 
angles in Fig. 4(b)) at the phase A-C and C-D (paramagnetic 
phase) phase boundaries. These experimental observations 
clearly suggest that Ba3NiNb209 evolves through successive 
magnetic phase transitions from the A phase to uud (B) phase 
and then to oblique (C) phase with increasing magnetic field. 
Whether a magnetization plateau should exist in the B phase is 
still an open question. Experiments performed on single crys- 
tal samples with in-plane magnetic fields at lower tempera- 
tures are expected to reveal the existence of the magnetization 
plateau more clearly. 

For TLAFs with classical spins, thermal fluctuations stabi- 
lize the uud phase at finite temperatures showing a plateau- 
like magnetization, hence the magnetic field range of the 
uud phase decreases with decreasing temperature. In 
Ba3NiNb209, the uud phase has a range of about 8 T even at 
0.5 K (<C Tni) and the width of the uud phase increases with 
decreasing temperature. This suggests that the uud phase has 
a quantum mechanical origin due to the small spin 5=1 for 
Ni^+. This makes Ba3NiNb209 a rare S= \ TLAF showing a 
uud phase driven by quantum fluctuations. Another notewor- 
thy feature in the phase diagram of Ba3NiNb209 is that within 
a narrow temperature window, 3.5 K < T < 4.5 K, with in- 
creasing field the spin configuration goes from the A phase 
to the C phase directly without passing through the B phase 
(see online supplemental material for details). This is signifi- 
cantly different from the other TLAF systems with uud phase 
mentioned in the introduction, either for classical or quantum 
spins. For those systems, in order to reach the C phase from 
the A phase with increasing field, it is always necessary to 
pass through the B phase. To our knowledge, so far no theory 
predicts this feature and we leave the explanation for it as an 
open question for future investigations. 



More interestingly, as the polarization measurements show, 
all these three magnetic phases have a FE ground state. Stud- 
ies of multiferroic properties of TLAFs have so far been lim- 
ited to systems with classical spins. Two celebrated examples 
are (i) ACr02 (A = Cu and Ag) with trigonal R-3m structure 
and (ii) RbFe(M04)2 with trigonal P-3m structure. It has been 
proposed that the ferroelectricity in ACr02 is driven by the 
helical spin-spiral order of the Cr'^^ {S = 'ill) spins among 
the stacked Cr-layers 11511 . The zero field magnetic structure 
of RbFe(M04)2 has AFM spins between the nearest neighbor 
layers and a 120 degree AFM ordering within the triangular 
layers. A phenomenological model based on symmetry argu- 
ments (the broken inversion symmetry and the absence of a 
mirror plane perpendicular to the c axis) was proposed to ex- 
plain the ferroelectricity in the A phase for RbFe(M04)2 ll4ll . 
A similar situation could lead to the ferroelectricity in the A 
phase for Ba3NiNb209, since the two Nb layers between the 
adjacent Ni layers lead to the absence of a mirror plane per- 
pendicular to the c axis [Fig. 1(e)]. Although a smaller po- 
larization, P = 1.1 i^lCIvp?, was observed due to the polycrys- 
talline nature of Ba3NiNb209, its value is still comparable to 
that of RbFe(M04)2 obtained with single crystal samples. 

A difference between RbFe(M04)2 and Ba3NiNb209 is 
the electric ground state at higher magnetic fields. In 
RbFe(M04)2, the ferroelectricity disappears in the B phase 
{uud phase), which was interpreted as the inversion symme- 
try in the triangular uud magnetic structure being reinstated. 
In contrast, the ferroelectricity persists in the uud and the 
oblique phases of Ba3NiNb209. The inset of Fig. 3(b) shows 
the magnetic field dependence of the polarization at 2.3 K and 
3.8 K obtained from the pyroelectric current measurement un- 
der ZFC. The polarization increases initially up to = 2 T fol- 
lowed by a monotonous decrease with the increasing magnetic 
field up to 9 T. Then, there is a change of curvature around 9 T, 
above which the polarization is larger than the projected value 
extrapolated from the linear behavior of the lower magnetic 
fields. This behavior suggests that the mechanism for the FE 
state in B and C phases is different from that for the A phase 
under lower magnetic fields. 

One possible explanation for the FE phase at high magnetic 
fields in Ba3NiNb209 is that the inversion symmetry is bro- 
ken by deforming the triangular lattice via magnetostriction. 
In the triangular spin structure, magnetostriction driven by 
spin correlation (Si • Sj) cannot break the inversion symmetry, 
but the anisotropic magnetostriction driven by external mag- 
netic fields can. In the distorted triangular lattice, the space 
inversion and the three fold rotation are not symmetric oper- 
ations any more, which can lead to finite polarization in the 
ab plane (see online supplemental material for details). Mag- 
netostriction often plays an important role in inducing mag- 
netic multiferroicity |23, l^l- The polarization values close to 
the B-C phase boundary at fixed T are too small to draw any 
conclusions, but the continuous polarization increase with de- 
creasing temperature while crossing the B-C boundary at fixed 
magnetic fields under the ZFC condition [Fig. 3(d)] suggests 
that the FE order parameter is same in the B and C phases. 
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From these observations, we speculate that the magnetostric- 
tion starts to play an important role in inducing ferroelectricity 
at high magnetic fields (B and C phases), while the spin chi- 
rality in the A phase is the driving force of the multiferroicity 
at low magnetic fields. 

In summary, although more work is needed to confirm the 
magnetic structure of the A, B and C phases, the detailed 
magnetic and electric studies on Ba3NiNb209 reveal that the 
system (i) is a new spin-1 TLAF system showing two mag- 
netic phase transitions bracketing an intermediate up-up-down 
phase driven by quantum fluctuations; (ii) is a multiferroic 
system with FE ground states stabilized at all three magnetic 
phases. These findings show that Ba3NiNb20g is a rare TLAF 
with strong couplings between the successive magnetic phase 
transitions and the ferroelectricity. One step further, we also 
confirmed the prediction that a broad range of trigonal mate- 
rials with the 120 degree spin structure can be multiferroic. 

This work was supported by NSF-DMR-06541 18 and the 
State of Florida. P. S. is supported by the DOE under Grant 
No. DE-FG02-98ER45707. 
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